Search for the photo-excitation of exotic mesons in the tt+tt+tt system 
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A search for exotic mesons in the tx^-k'^-k~ system photoproduced by the charge exchange reaction 
7P — > 7r"'"7r+7r~(n) was carried out by the CLAS collaboration at Jefferson Lab. A tagged-photon 
beam with energies in the 4.8 to 5.4 GeV range, produced through bremsstrahlung from a 5.744 
GeV electron beam, was incident on a liquid-hydrogen target. A Partial Wave Analysis (PWA) 
was performed on a sample of 83,000 events, the highest such statistics to date in this reaction at 
these energies. The main objective of this study was to look for the photoproduction of an exotic 
J^*^ = 1 ^ resonant state in the 1 to 2 GeV mass range. Our PWA analysis, based on the isobar 
model, shows production of the a2(1320) and the 7r2(1670) mesons, but no evidence for the ai(1260), 
nor the 7ri(1600) exotic state at the expected levels. An upper limit of 13.5 nb is determined for the 
exotic 7ri(1600) cross section, less than 2% of the a2(1320) production. 

PACS numbers: 13.60.Rj, 13.60.-r, 14.60.jn 



The self interacting nature of the gluon within Quan- 
tum Chromodynamics (QCD) allows for hybrid resonant 
states with a {qqg) configuration, where the gluonic de- 
gree of freedom gives rise to a spectrum of additional 
states outside the constituent quark model (CQM). The 
observation of gluonic-hybrid hadrons, with an explicit 
excitation of this gluonic degree of freedom, will be an 
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important test of the predicting power of QCD at inter- 
mediate energies. A signature of states beyond the CQM 
would be the existence of mesons with quantum num- 
bers that cannot be attained by [qq) mesons (so-called 
"exotics"). In particular, the lowest lying {qqg) state is 
predicted to have J^^ = 1^+ [H, d, [1,11] ) and mass near 
1.9 GeV. A more recent calculation on the lattice using 
lighter quark masses predicts a mass at 1.74 GeV Q. 

Here we report on a search for exotic mesons decaying 
to three charged pions. This channel was chosen for its 
simplicity, since only a few decay channels are available 
to this final state. In addition, although the dominant 
decay mode of the lightest J^'~' = 1 ^ state is predicted 
to be into an S-wave and P-wave meson, such as the 
&i(1235) TT or /i(1285) tf the three-pion final state in 
the p decay channel could be non-negligible 0, H, H] ■ 

There is also evidence for an exotic 7ri(1600) state 
in the reaction 7r~p TT~^Tr~^Tr~p at 18 GeV by the 
Brookhaven E852 experiment [H, [ll| . Furthermore, the 
existence of the 7ri(1600) has been confirmed by the same 
experiment in other channels such as ryV T^. The VES 
collaboration also reported an exotic signal in the ryV^ 
channel in the same mass region [l3| . However, a more 
recent analysis of a higher statistics sample from E852 
Stt data, claims to find no evidence for 7ri(1600) produc- 
tion [l3|- AH these results were obtained by pion beam 
experiments. 
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It has long been anticipated that photoproduction ma' 
be a better production mechanism for exotic mesons 
i, [H, [H, [13. If the 7ri(1600) state couples to pir, 
then this state should also be produced with a photon 
beam through tt exchange via Vector-Meson-Dominance 
(VMD) [l8|. However, the use of photon beams as probes 
for exotic meson production has not been fully explored 
so far, and the existing data on multi-particle final states 
are sparse. Photoproduction of mesons in multi-pion final 
states has been reported in three previous experiments: 
a SLAC 40-in. hydrogen bubble chamber experiment [l^ 
at 4.3 GeV, a CERN hydrogen experiment [2^ between 
25 and 70 GeV, and a SLAC 1-meter hybrid bubble cham- 
ber experiment [2l|, |22| around 19 GeV. All these ex- 
periments lacked the statistics required for a full Partial 
Wave Analysis. For example, Condo et al. 22] showed 
that the three-pion spectrum in the low-mass region (be- 
low 1.5 GeV) is dominated by a2(1320) production, and 
found no clear evidence for ai(1260) production. In the 
high-mass region (1.5 to 2 GeV), based on an angular dis- 
tribution analysis, this same group claimed evidence for 
a narrow state at 1.775 GeV with possible J^*^ = 1 ^, 
2 '", or 3"'""'" assignments. From the analysis of four pion 
events, Ref. [2^ reported two peaks, one at the mass of 
the 02(1320) and another at around 1.75 GeV, in the pn 
state recoiling off the proton and the remaining pion. 

The present experiment was performed at Jefferson 
Laboratory during the 2001 running period of CLAS. 
Details on the design and operation of CLAS and its 
components may be found in Ref. p3j and references 
within. The experiment ran with a fixed electron beam 
energy of 5.744 GeV. A tagged-photon beam, with a flux 
of 5 X 10^ photons/s produced via bremsstrahlung from 
a 3xl0~^ radiation lengths radiator, was incident on a 
liquid-hydrogen target contained in a cylindrical cell 18 
cm in length. The running conditions were optimized 
for meson production recoiling off a neutron: the CLAS 
torus magnet field was reduced to half its maximum cur- 
rent, and the target pulled back 1 m from the nominal 
center of CLAS position. Only events with tracks com- 
ing from the target material that were in time within 1 
ns with a beam photon with energy between 4.8 and 5.4 
GeV were selected. 

The three pions in the reaction 77? — s- 7r+7r+7r~n were 
identified by the time-of-flight detector, and the neutron 
was identified through missing mass (Fig. [1]). The re- 
constructed mass of the neutron was 942.2 MeV with 
a width of 25.1 MeV. As shown in Fig. [^h, ^J, the vari- 
ous nTT invariant mass spectra show clear production of 
baryon resonances recoiling off a two-pion system, consis- 
tent with the process shown in Fig.[3f). The titt^ invariant 
mass (Fig.[5]i) shows a peak at the mass of the A(1232), 
while the n7r+ effective mass distribution (Fig. ^) shows 
enhancements around the mass of N*(1520/1535) and 
N*(1650/1675/1680), with a smaller contribution from 
A(1232) production. 

The observed squared four-momentum transfer distri- 
bution, — i, between the incoming photon beam and the 
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FIG. 1: Missing mass of the three-pion system. The distri- 
bution was fit to a gaussian plus a linear background (solid 
curve). The inset shows a wider scale. The region between the 
vertical lines at 0.884 and 0.992 GeV identified the neutron. 
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FIG. 2: Mass spectra, (a) mx~ invariant mass, (b) nvr^ in- 
variant mass (the left arrow indicates N*(1520/1535) reso- 
nances and the right arrow N*( 1650/ 1675/ 1680) resonances), 
(c) TT^TT^ invariant mass, and (d) -k^-k^tx" invariant mass. 
The hatched areas are the distributions of remaining events 
after the momentum transfer and tt"*" laboratory angle selec- 
tions. The histograms in (b) and (c) are filled twice for each 
event, once for each tt^. 



produced three-pion system follows an exponential of the 
form f{t) = ae~''l*l. This form is consistent with periph- 
eral production of the three-pion system recoiling off the 
neutron and is consistent with the production mechanism 
shown in Fig. For the exponential slope constant, 6, 
we obtained a fitted value of 3.02 GeV^^. To enrich 
the sample of peripheral events, we required \t\ < 0.4 
GeV^. To further enhance the mesonic sample relative to 
the baryonic background, only events with forward-going 
TT+'s (0^06(71"*") < 30°) were selected. This requirement, 
as determined from simulation, is necessary to remove 
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most of the baryonic background without significantly 
affecting the mesonic sample. After all the former se- 
lections were applied, about 83,000 events remained as 
input to the Partial Wave Analysis (PWA). 
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FIG. 3: (a) Peripheral meson production, (b) baryon reso- 
nance background processes. 

The TT+TT^ invariant mass distribution. Fig. ^jp, shows 
enhancements at the mass of the p(770) and /2(1270). 
The Stt invariant mass spectrum. Fig. [2}i, shows a clear 
peak around the mass of the 02(1320) meson. It also 
shows a broad enhancement in the 1600 - 1700 MeV mass 
region, which could be due to the 7r2(1670) or 7ri(1600), 
but it could also be the result of acceptance and phase 
space limitations. 

A Partial Wave Analysis, based on the isobar model, 
was performed to determine the spin, parity, and charge- 
conjugation (J^^) of the three-pion system. The three- 
body decay of the meson was described by a sequence 
of two-body intermediate decays (Fig. States with 
definite J^^ of the mesonic system were combined as 
eigenstates of reflection through the production plane 
(reflectivity) according to the formalism described in 
Refs. ^2^, i25|, i26|, i27,] . The expected number of events per 
mass bin were written as a sum of production and decay 
amplitudes. To obtain the observed number of events 
per bin, the acceptance of the detector was determined 
by normalizing to the Monte Carlo sample. The model 
used the same t distribution as determined from the data, 
and the events were simulated in a GEANT-based model 
of the CLAS detector. A maximum likelihood fit to each 
mass bin was done with a set of input partial waves. 
Many different sets of waves were tried during the anal- 
ysis. The meson spin-density matrix was approximated 
by a rank 2 matrix due to limited statistics and the large 
number of parameters required for the fit. The rank of 
the spin-density matrix represents the number of inde- 
pendent spin components in the initial and final states 
necessary to describe the interaction. Parity conserva- 
tion reduces the number of independent components to 
four. 

No baryon resonance partial waves were used and an 
isotropic non-interfering background wave was included, 
which is appropriate, since the events were dominated 
by meson production after the kinematic selection cuts. 
Below 1.36 GeV, no /27r wave was included as its mass 
is lower than the nominal threshold. Contributions from 
TO^ = 0+,0^ waves were found to be small, and these 
waves are therefore not included in the final wave set. 
The final list of waves that provided the best stable fit 
is shown in Table [J The quality of the PWA fit was 



checked by comparing various mass and angular distri- 
butions from the observed data and the predictions of 
the PWA results, i.e., weighting the Monte Carlo events 
according to the PWA fit. Qualitatively good agreement 
was observed between these two sets [28l] . 
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TABLE I: Set of partial waves used in the PWA. J, P, C refers 
to the Spin, Parity, Charge-conjugation of the three-pion sys- 
tem, m, e are the spin-projection and reflectivity quantum 
numbers. L is the orbital angular momentum of the decay. 
Isobar is the intermediate meson that decays to two pions. 
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FIG. 4: PWA results: Combined intensities of waves included 
in the flt. The intensities shown are the sum of the intensities 
from both ranks of the spin-density matrix: (a) Total inten- 
sity, (b) Background intensity, (c) 2'^'^{ptt)d intensity, and 
(d) 2""*" total intensity. The dashed vertical line in (c) and 
(d)) indicates the border of the two fitting regions. 
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Results from the PWA are shown in Figs. [H and [S] 
The two most prominent waves are J^*^ = 2++ and 2 ^ 
{F\g.\^,d). The 02(1320) is observed in the {ptt)d wave 
at the expected mass and width. For the 02 cross section, 
we obtained a value of 0.81 ± 0.25 The uncertainty 
includes both detector systematics and variation of PWA 
results estimated by using different wave sets and starting 
parameters. The result is in rough agreement with the 
previous measurement of 1.14 ± 0.57 jjh obtained at 5.25 
GeV by Ballam [H and 1.71 ± 0.86 ^lh obtained at 4.7 to 
5.8 GeV by Eisenberg 

The 7r2(1670) is observed clearly in the {f 21^)3 de- 
cay mode (Fig. [5t»). In the remaining 1^2 decay modes: 
(/97r)p_i?(Fig.[SJi), and {f2T^)D, the signal looks broad and 
distorted in the low-mass region. This is most likely due 
to the incomplete set of waves in the fit, and is indica- 
tive of leakage from the 02. Similar observations were 
reported by the E852 collaboration [O, [l3]- In addition, 
due to the two different wave sets used in the two fit- 
ting regions (Table some instability is observed near 
the border of 1.36 GeV (Fig.®, d). This should be ex- 
pected as a result of incomplete waves sets and can only 
be resolved by future higher statistics experiments [3l|. 

It is important to note that we see no clear evidence 
for the ai(1260) in the possible decay mode of the {pTr)s 
wave ( Fig . \5jp). This observation agrees with Condo 
et al. [2l|. The event enhancement observed around 
1.3 GeV in the non-interfering, flat background waves 
(Fig. is likely leakage coming from the 02 and may 
demonstrate that the rank of the spin density matrix is 
larger than 2. 

We do not observe resonant structure in the exotic 
1 ^(/57r)p partial wave (Fig. [SJi). To determine if the 
1 ^ wave was necessary to better fit our data, we com- 
pared PWA fits with and without the 1 ^ wave using 
the Likelihood Ratio (LR) test [H, [s^l . Using LR statis- 
tics, we find that the PWA set of waves including the 
1 ^ wave fits the data significantly better than a model 
without the 1 wave. While no clear resonant structure 
was observed in the 1 intensity distribution (Fig. [5]i), 
this distribution was used to estimate an upper limit to 
the 7ri(1600) cross section using the method of O. Hc- 
lene [34|. Using the mass of 1597 MeV and the width of 
340 MeV as measured by Ref. [l2| , we estimated an upper 
limit for the tti (1600) of 13.5 nb at a 95% confidence level, 
less than 2% of the 02(1320). Therefore, our results do 
not agree with the predicted strengths for pho toproduc- 
tion of a 1 gluonic hybrid meson [gI. [sI. [l5l.]T7l? ]. Based 



on Ref. [1], the 7ri(1600) is expected to be produced with 
a strength near 10% of the 02(1320). Ref. |T3| predicts a 
factor of 5-10 larger ratio of exotic meson to 02 in photo- 
production than hadroproduction. This would imply the 
7ri(1600) cross section to be on the order of 50% of 02, 
more than 25 times higher than what we observed. It is 
possible that the 7ri(1600) reported by Ref. [l3| is not a 
gluonic hybrid meson, but rather of other nature. Alter- 
natively, the calculated photoproduction cross section of 
gluonic exotic mesons might be overestimated. 
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FIG. 5: Partial wave decompositions: (a) 2 {p-k)p^f (the 
dashed vertical line indicates the border of the two fitting 
regions), (b) 2- + (/27r)s, (c) l++(p7r)s, and (d) 1-+{ptt)p ( 
the vertical and horizontal arrows indicate the mass and width 
of the 7ri(1600) as reported by Ref. [iC 
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